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We present results for the imaginary parts of the I — and 2 K — > mr decay amplitudes, Irm4o,2, and the ratio 
of CP violation parameters, e'/e. Our calculation is done in the quenched approximation at a -1 ~ 2 GeV, lattice 
size 16 3 x 32, using domain wall fermions with L 3 = 16. We study the three flavor case (charm is not an active 
flavor) and find e'/e i s small and slightly negative. 



1. Introduction 

In this talk we present results for the CP vi- 
olating imaginary parts of the isospin zero and 
two K — > 7T7T decay amplitudes, and hence the 
theoretical value of e' . We build on the talks by 
Mawhinney and Cristian where details of the ba- 
sic method and simulations can be found. Cris- 
tian presented results for the CP conserving real 
parts of the amplitudes, and it was found that 
the theoretical calculation of the ratio of the 
real parts was in good agreement with the long- 
standing AI =1/2 rule. Here we find the imag- 
inary parts, together with the real parts, yield 
a value of e'/e = (-4.0 ± 2.3) x 10~ 4 (statisti- 
cal errors only) that is small and slightly nega- 
tive, in contrast to the recently determined val- 
ues, - (15-20) x 10" 4 , from the NA48 and KTeV 
experiments at CERN and Fermilab. We refer 
the reader to |l| for complete details of our cal- 
culation. The CP-PACS collaboration has also 
presented a very similar calculation at this meet- 
ing @ and finds results for e'/e that are compati- 
ble with ours. 

Defining decay amplitudes as 

A(K° -> 7T7r(/)) = Aie iSl , 

A(K° inr(I)) = ~A}e tSl , (1) 
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where / = 0, 2 is the isospin of the final state 
pions and Si the corresponding s-wave phase shift, 
it can be shown that 

e' _ i e-^o-^) r c a 2 /ImA 2 ImA \ 
7 _ V^e Re A \Re A 2 ~ ReA / ' 

where we have neglected isospin breaking effects 
(these are expected to decrease the value of e'). 
The factors outside of the parentheses are well 
known experimentally, so we quote values of e'/e 
using these known values. Thus, the calculation 
boils down to determining the difference of ratios 
Pj = ImAr/ReAr. 

Initio an d ImA 2 are given by Im(e~ l(5l x 
{■Kir{I)\-m { ^ s=1) \K )) where the effective weak 
Hamiltonian, TL^ AS=1 \ is a linear combination of 
four-quark operators with Wilson coefficients de- 
termined by the Standard Model. The imaginary 
part of 7i( AS=1 ) arises from the CP odd phase of 
the CKM mixing matrix. The non-perturbative 
calculation of the matrix elements of the four- 
quark operators forms the fundamental core of 
our calculation (see the talk by Mawhinney). We 
use chiral perturbation theory to extract the low 
energy constants of 7i( AS,=1 ) from calculations of 
unphysical K — > 7r and K — > matrix elements. 
This procedure directly yields the corresponding 
K — > mr amplitudes to lowest order in chiral per- 
turbation theory || (O(p , p 2 ) depending on the 
particular operator). Some of the 0(p A ) correc- 
tions to these relations have been calculated , 
which we include in our analysis. However, a com- 
plete next-to-leading order calculation using only 
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the reduced matrix elements cannot be done since 
the relevant 0(p 4 ) counter- terms in the AS = 1 
chiral Lagrangian differ in the two cases, ft must 
therefore be argued that these higher order effects 
are small in the final answer. With the low en- 
ergy constants and Wilson coefficients in hand, it 
is straightforward to calculate the physical decay 
amplitudes. 

2. Results 

In Fig. |f] we show ImAa as a function of a 
fictitious mass parameter £ {m 2 K , m 2 . — ► £ x 
(mfc, n%„)) which we have introduced in order to 
study the behavior of the "physical amplitudes" 
from the chiral limit (£ = 0) to the actual phys- 
ical point (£ = 1). In particular, our calculation 
should be increasingly accurate as £ — > 0. We 
show results using lowest order chiral perturba- 
tion theory as well as those with one-loop correc- 
tions. The amplitude is dominated by the QCD 
penguin operator Qq which vanishes in the chi- 
ral limit, ft does not vanish at £ = due to the 
small contribution of the electroweak penguin op- 
erators Qr t 8. We note that the one- loop chiral log 
corrections are significant (roughly 50%). 

Similarly, in Fig. || we show ImA2. Since Imyl2 
is dominated by Q% whose matrix element is a 
constant at lowest order, it is also roughly con- 
stant over the whole range, and the one-loop cor- 
rections are not as large as for the 1 = channel. 
It is interesting to note that at the physical point 
the imaginary parts of the amplitudes also exhibit 
a "A/ = 1/2" rule similar to the real parts. 

The residual renormalization scale dependence 
in the physical amplitudes is shown in Fig. ^| and 
is most apparent for the 1 = 2 case where the 
statistical errors are smaller. In both cases, the 
largest deviation is at the lower end. Had the 
Wilson coefficients been calculated to all orders 
in perturbation theory, had our matrix elements 
been evaluated in full QCD, and had we included 
subtractions of all operators with lower dimen- 
sion (not just power divergent ones) this scale de- 
pendence would have canceled exactly in 7i AS=1 
since its matrix elements are physical observables. 
As fj, decreases non-perturbative effects in the op- 
erators may arise which are not compensated in 
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Figure 1. The imaginary part of Aq. One loop 
chiral log corrections are large at the physical 
point, £ = 1. The errors are statistical only. 



the Wilson coefficients. On the other-hand, if /i 
becomes too large, lattice artifacts appear in the 
renormalized operators. Thus we quote values for 
physical quantities at fi = 2.13 GeV (of course, 
the scale is essentially determined by the inverse 
lattice spacing, a -1 = 1.922 GeV in this case). 

Our results for e'/e are shown in Fig. ^. 
The value at the physical point is small and 
slightly negative whereas the experimental value 
lies somewhere between 15 and 20 parts in 10~ 4 . 
Interestingly, the one-loop corrections that were 
significant in the individual amplitudes do not 
play any role here. Just away from the chiral 
limit, both the numerator and denominator of Pq 
are dominated by (8,1) operators Qe and Q2, re- 
spectively, so the log correction (which is the same 
in both cases) cancels. Since the leading behav- 
ior of both is £,{m 2 K — m 2 ), Pq behaves essentially 
like a constant (see Fig. Q). P2 is dominated by 
Qs in the numerator and Q2 in the denominator. 
The corresponding one-loop corrections are not 
the same; however, they are similar in magnitude, 
so the net effect is not large. Since Re^ is dom- 
inated by (27,1) operators, it (almost) vanishes 
as £ — ► whereas we saw Im^ was roughly con- 
stant. So, as £ — > 0, e'/e becomes large and neg- 
ative (note that at £ = 0, Pq = P2 and e'/e = 0). 
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Figure 2. Same as Fig. [l] except for Ieou4.2. 



3. Summary and Outlook 

We have completed a first calculation of e'/e 
in the Standard Model, albeit with significant 
approximations including quenching, use of per- 
turbation theory below m c harm, and incomplete 
next-to-leading order calculations in chiral per- 
turbation theory. Such approximations could 
conspire to yield the current result which is in 
disagreement with experiment, even though it is 
difficult to point to a single approximation that 
may be responsible. In any case, the above ap- 
proximations are not insurmountable and will be 
addressed in future calculations. In particular, 
the results reported here are the first part of a 
two-part calculation. The second part, not yet 
finished, uses an active charm quark and the four- 
quark effective theory. The road ahead will be 
interesting, especially if we must turn away from 
the Standard Model. 
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Figure 3. Residual rcnormalization scale depen- 
dence in the physical amplitudes (units are GeV). 
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Figure 4. The ratio of CP violation parameters 
e'/e. Totals are shown for lowest order (upper) 
and one-loop (lower) chiral perturbation theory. 
The 1 = and 2 contributions are shown for the 
one-loop case only. The experimental value is be- 
tween 15 and 20. 



